Screening in anisotropic superfluids and the superfluid density in underdoped cuprates 
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We examine the nature of the collective excitations in a strongly anisotropic system of bosons 
interacting via Coulomb interaction. Such a system has often been used in the past to model 
the effects of quantum and classical phase fluctuations on the superfluid density of underdoped 
cuprates. Depending on the anisotropy and the effective strength of the interaction we find four 
different regimes for the temperature dependence of the superfluid density. Coulomb interaction in 
underdoped cuprates is argued to be effectively short-ranged, and less then unity in appropriately 
defined units. 



c 

o 
o 

I 

Or 

■4-* 



I 

-a 
c 

o 
o 



> 

VO 
(N 

O 

in 
o 

> 



i 

C 

O 
o 



Temperature dependence of the superfluid density 
p(T) is one of the key properties of a superfluid. This is 
particularly true in high-temperature superconductors, 
where the linear low-temperature behaviour of p(T) at 
optimal doping served as early evidence for the d-wave 
symmetry of the superconducting order parameter 0. 
Another intriguing aspect of p(T) is its evolution with 
underdoping, the effect of which may be expected to in- 
crease the importance of fluctuations. Indeed, it has been 
argued that for a low value of p(0) as is observed in the 
underdoped regime the further reduction of p{T) with 
temperature should be primarily due to classical phase 
fluctuations 0, 0- This argument has been criticized 
for its neglect of Coulomb interaction, which is expected 
to strongly suppress classical phase fluctuations below 
the plasmon energy and thus reinstate quasiparticles as 
the source of the linear temperature dependence of p(T) 
at low T 0. The competition between the quasiparti- 
cles and phase fluctuations for the form of p(T) becomes 
quite explicit in the so-called Ioffe-Larkin rule 5] 
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P-\T) = p~ p \T) + pjl u jT), 



(1) 



where p qp {T) is the standard BCS quasiparticle contri- 
bution in a d-wave state, and pfiuctiT) is the (model- 
dependent) fluctuation component. This transparent re- 
sult was first derived within the context of effective gauge 
theories of the t- J model 6] , but it may be expected to 
apply more generally to strongly fluctuating quasi-two- 
dimensional superconductors with p(0) <C p qp (0) pfl- 

Cuprates, however, are strongly anisotropic materials, 
with the anisotropy between the transport properties in 
the ab-plane and the c-axis increasing with underdoping. 
It is well known that under these conditions the plas- 
mon dispersion becomes very anisotropic as well, and the 
large plasmon energy gap gets replaced by a much lower 
one proportional to the interlayer coupling 0, 0. One 
may therefore expect Coulomb interaction in such highly 
anisotropic superconductors to become less efficient in 
gapping the phase mode. In this paper we study in 
greater detail the combined effect of large anisotropy and 
Coulomb interactions on p(T). We model the fluctuation 
component in Eq. ^ by a layered system of bosons in- 
teracting via Coulomb interactions and with the density 
proportional to doping. Such an effective theory arises 
naturally in several theories of underdoped cuprates , 
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FIG. 1: Four regimes for the temperature dependence of the 
superfluid density p(T) in a layered bosonic system with weak 
Coulomb interactions: I) quasi-two-dimensional (2D) regime 
with weak and effectively short-range (screened) interaction 
for K(Q)/t 3> 1/A, II) quasi-2D regime with weak long-range 
interaction for 1 <C K(0)/t <C 1/A , III) three-dimensional 
(3D) regime with weak long-range interaction for A 3 <C 
K(0)/t <C 1, and IV) 3D regime with strong long-range inter- 
action for K(0)/t < A 3 . K(0) = h 2 p(0)/m, and t is the inter- 
layer Josephson coupling. A = (2ne 2 /(ed))/(h 2 /(md 2 )) <C 1 
is the dimensionless strength of the Coulomb interaction, and 
d is the inter-layer separation. 



0, [13, d, and provides a rather general representa- 
tion of a charged layered superfluid. We begin with the 
simple case of a two-layer system and demonstrate that, 
when the layers are Josephson-decoupled, there exists a 
linearly dispersing phase mode at low wavevectors. This 
is essentially a consequence of the perfect screening of in- 
teractions in one layer by fluctuations in the other. Gen- 
eralizing to a system with an infinite number of layers, 
for a weak Coulomb interaction we find that there are 
four discernible regimes for the temperature dependence 
of the superfluid density, controlled by the ratio between 
the boson density and the Josephson coupling (Fig. 
An estimate of the relevant parameters places the under- 
doped YBa 2 Cu 3 6+a; (YBCO) with T c > 5K firmly in 
the regime I, where the Coulomb interaction acts as an 
effective short-range interaction; the long-range nature of 
the Coulomb interaction in this regime is irrelevant ex- 
cept at extremely low temperatures. We also find that 
the value of the dimensionless interaction strength A in 
YBCO is A < 1, and argue that it is possible that A <C 1. 

Let us define the quantum mechanical action for a lay- 
ered system of interacting bosons as S = f dr J d 2 xC 
with 
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C = 5>*(x,r)(d r 



2m 



(2) 



J d 2 x , \^ i (x,r)\ 2 V ij (x-x')\^(x',rf 



where iV is the number of two-dimensional (2D) layers, 
i is a weak Josephson coupling between the layers, p = 
1/T, and we set ks — h — 1 throughout, unless otherwise 
noted. The Coulomb interaction is 

V l3 {x-x') = e 2 /eyJ\x-x'\ 2 + \i-j\ 2 d 2 ), (3) 

d being the separation between the layers and e the static 
background dielectric constant. We assume the presence 
of a neutralizing background of density po equal to the 
average areal density of bosons at the chemical potential 

The nature of the excitations in the above superfiuid 
system is most explicit in a system of only two layers, 
which differs already from a single layer in an impor- 
tant way. Let us introduce first the usual density-phase 
variables [Tg| as $,(x,r) = \f po + Ui(x, r)e iei ^' r) and 
expand the Lagrangian to the quadratic order in IT and 
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-t Po (e 1 -6 2 ) 2 + -^(n 1 -u 2 ) 2 . 



We assumed here the usual periodic boundary conditions 
in imaginary time: Hi(x,/3) — n,(a?, 0) and 9i{x,0) = 
9i(x, 0) + 2Trrii(x), with ni(x) integer. By rotating the 
fields as 9± = {9 1 ±9 2 )/V2, U± = (LTi ±U 2 )/V2, we can 
decouple the above Lagrangian as C = £+ + with 

c + = £(v^) 2 + ^ + n + ). + + ^( 5 ) 

+i J d 2 x'n + (x,r)K f (f-x')n+(x',r), 

(vn_) 2 



£- = ^(V9_)^ + m_( 
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(6) 



+ 2^ 2 + -^n 2 _, 

where V±(x) = Vu(x) ± Vi2(:r). The Gaussian integra- 
tion over II± yields two branches of excitations with the 
energies 
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FIG. 2: The two branches u>± of the excitation spectrum of the 
two-layer system, with t = and weak Coulomb interaction. 
The lowest mode crosses over from linear behaviour at low k 
to k 2 behaviour. The plasmon starts out as \fk before crossing 
to k 2 . 



where 
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2ire 2 
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The branch u + describes the usual two-dimensional plas- 
mon, uj + w yj Aire 2 pok / em at low momenta. The two 
layers oscillate in phase and, as a consequence, lu + is in- 
dependent of the Josephson coupling. The canonically 
conjugate variable to 9 + is the sum of two densities, and 
therefore the Coulomb interaction affects the energy of 
this mode, as is usual in systems with long-range inter- 
actions. In contrast, the conjugate variable to 9- is the 
difference between the two densities which can oscillate 
without any cost in Coulomb energy. As a result, for 
t = 0, uj- w (2npoe 2 d/em) 1 / 2 k at low momenta, and the 
dispersion of the lower branch is the same as if the system 
had only a short-range interaction of strength ~ 2ire 2 d/e. 
With t = the layers cannot exchange particles and the 
density in each layer therefore may oscillate so as to per- 
fectly screen the Coulomb interaction in the other. The 
oscillations, however, are then out of phase, and conse- 
quently when t =/= this mode becomes gapped, with 
cij_ sa y / 2ntpoe 2 d/e at low momenta. We will refer here- 
after to this energy as the Josephson gap. 

The remarkable feature of the above result is that in 
a system with negligible Josephson interaction between 
the layers, the Coulomb interaction becomes effectively 
short-ranged as far as the low-energy excitation spec- 
trum is concerned. More precisely, when t — 0, w_ is 
linear and deviates from lo+ significantly for k -C 1/d. 
For a large separation between the layers and for 1/d <C 
k <C (87rm/9 e 2 /e) 1 / 3 , w_ w uj + ~ vfe. Finally, for 
(87TTOp e 2 /e) 1/3 < k, lu- « uj+ « k 2 /2m. If, on the 
other hand, the interaction is weak or if the layers are 
brought close together so that 1/d ^> (8irmpoe 2 /e) 1 ^ 3 , 
w_ ~ k for k -C y^8irmpoe 2 d/e, and w_ ss k 2 /2m other- 
wise, without the intermediate region oj_ ~ vfc. In this 
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regime w_ becomes identical to the phonon spectrum of 
the weakly interacting Bose gas. This is illustrated in 
Fig. With a finite Josephs on c oupling w_ approaches 
the Josephson gap for k <C y/Amt. If the Josephson cou- 
pling were strong, of course, a->_ ^> lo+, and the plasmon 
would resume its place as the low-energy mode of the 
system. 

In a system with N layers, for t = there are TV — 1 
modes with linear dispersion and only a single plasmon. 
This is easily established by considering the interaction 
matrix V tJ {k). When k -> 0, V tJ {k) = (2ne 2 /ek){l + 
0(kd)). So in the limit k -> 0, 



V(k) 



2ne< 



(1,1, ...1) T 0(1,1,....!), 



(10) 



and the interaction matrix has one eigenvector with the 
eigenvalue 2irNe 2 /efc, and N— 1 degenerate eigenvectors 
with zero eigenvalue. The former eigenvector is the to- 
tal density which is canonically conjugate to the sum of 
the phases and describes the plasmon. The latter N — 1 
modes, being orthogonal to the plasmon, are electrically 
neutral and consequently cross from linear dispersion at 
low momenta to the Josephson gap at k = 0. 

The existence of linearly dispersing modes below the 
usual plasmon modifies the behavior of the superfluid 
density at low temperatures. To be specific, we focus on 
the system with infinitely many layers which is relevant in 
the context of high-temperature superconductivity. Im- 
posing periodic boundary conditions in the direction or- 
thogonal to the layers the excitation spectrum becomes 

Lu 2 (k, k z ) = e(k, k z ){2 Po V(k, k z ) + e(k, k z )), (11) 

with e(k, k z ) = (k 2 /2m) + tsm 2 (k z d/2), and [3 

2ne 2 sinh(kd) 



V{k,k z ) 



ek cosh(fceJ) — cos(k z d) 



(12) 



For k z = one finds the usual three-dimensional plasmon 
at w 2 (0, 0) — lo 2 = Aire 2 po/ dmc, while when k z ^ and 
t = 0, for k <C 1/d, u}(k,k z ) — uj p k/k z . The latter modes 
become gapped when t ^ and u>(k — + 0, k z ^= 0) = 
ojp^ftmd 2 /2. 

The temperature dependence of the areal in-plane su- 
perfluid density in Landau's two-fluid model is given by 



P(T) 



P(0) 
_d_ 

2m 



(13) 
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* /d dk z h2 dn b {u{k,k z )) 
-n/d 2 ^ duj(k,k z ) 



where n^iuj) is the usual boson occupation number |14| . 
We will find it convenient to express the superfluid den- 
sity in units of energy by defining K(T) = K 2 p(T) / m; we 
have also restored the dimensionful Planck constant. The 
rescaled superfluid density K{T) may then be expressed 
entirely in terms of dimensionless quantities as 



K{T) = K(0) 



T 
8^ 



(14) 



ydy / dzsinh 



2XK(0) smh\J2yT 



1/2 v 



T 3 / 2 V217(cosh J 2yT - cos(ttz)) , 



where f(y,z) = y + {t / T) sin 2 (irz/ 2); X = X/T d are 
dimensionless, with Td — h 2 /(md 2 ) as the characteris- 
tic energy scale in the problem. The parameter A = 
2we 2 / (edTj) is the dimensionless measure of the Coulomb 
interaction's strength. 

Equation ^5] or ^] expected to be valid for A <C 1 
and not too close to the critical temperature, leads to 
four distinct regimes of temperature dependence of the 
superfluid density. Take t/K(Q) <C 1, as is relevant to 
the cuprates, and consider the function K(T) as K(0) 
is decreased at fixed t. We will assume this to crudely 
correspond to underdoping a high-temperature supercon- 
ductor, as we discuss shortly. 

I) For K(0)/t 3> 1, the system is quasi-2D. Assuming 
A <C 1, to the zeroth order in A the superfluid density in 
Eq. is easily seen to equal the Bose condensate in the 
layered non-interacting system , and 



1.386 + 0(t/T)), (15) 



over most of the temperature range. The deviations from 
Eq. 15 are most significant below the crossover temper- 
ature T x = XK(0), where AK(T) = K(0) - K(T) - T 3 , 
and within the critical region of width ~ XT C around T c 
with T c w 2irK(0)/\n(2nK(0)/t). Besides T u there ex- 
ists also a lower crossover temperature T[ of the order of 
the Josephson gap, where AK(T) becomes exponentially 
suppressed. The latter temperature is 
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Xt , 2ttK(0) 
In ■ , 



whereas 



2tt V K(0) 



Tl _ X_ j 2ttK(0) 

Y c ~ 2^ 11 t ' 



(16) 



(17) 



The three characteristic temperature scales for variations 
of K(T) will therefore satisfy the inequalities 



T{ < Ti < T c 



for 



K(0) 



«A« 



2tt 



(18) 



(19) 



For such an interval for A to exist we obviously need 
t/K(0) < 27r/ln(27rif(0)/t), which is comfortably satis- 
fied for K(0)/t > 1. When the inequality l(TH|l is satisfied 
the long-range nature of the Coulomb interaction is rel- 
evant only at very low temperatures and the interaction 



appears in p(T) as being effectively short-ranged, and 
weak. Furthermore, as the ratio K(0)/t is reduced the 
relative temperature range over which K(T) behaves as 
a power law, Ti/T c , decreases as well, albeit logarith- 
mically slowly. One may interpret this decrease as that 
the interaction is being slowly renormalized towards zero 
with the reduction of the boson density 1 1 ol | . The regime 
in which the inequalities i|18|) hold is marked as I in Fig. 
□ 

II) If K(0) is decreased further so that K(0)/t ~ 1/A, 
one finds T[ m T\ -C T c . The low-temperature behavior 
of the superfluid density corresponds to the long-range 
interaction. This regime could be named the 2D, weakly 
interacting, long-range regime, and is labeled II in Fig. 
Ill 

III) For K(0)/t ~ 1, K(T) to the zeroth order in A, has 
the form of the fully three-dimensional Bose condensate 
over most of the temperature range, and 



K{T) « K(0) - 1.306 



T 3/2 



(20) 



In this regime the exponential behavior of AK(T) sets 
in below 



X 1 / 2 / 

The effective short-range behaviour AK(T) 
other hand, would appear below 

T 2 = AX(0)«^(^)Va Tc . 

7T t 



rp4 



(21) 



the 



(22) 



Note that TijT c now decreases as a power of K(0), which 
is again equivalent to the infrared renormalization group 
flow of the short-range coupling constant in a weakly 
interacting system of 3D bosons 0. However, since 
T2 > T 2 for K(0)/t < 1/A, and 1/A > 1, by the time the 
system enters the 3D regime where K(Q)/t ~ 1 the tem- 
perature dependence of K{T) crosses over directly from 
exponential at low temperatures to ~ T 3 / 2 at higher tem- 
peratures. This is then the 3D, still weakly-interacting, 
long-range regime, labeled III in Fig. ^ 
IV) Eventually, by reducing the density further one enters 
the regime where the Josephson gap becomes comparable 
to T c . From Eq. 21 this occurs when 



K(0) _ A 3 

t ~ 7T 6 



(23) 



This is the strongly-interacting regime (labeled IV in Fig. 
^) in which the superfluid density varies exponentially 
over the scale of ~ T c . In this regime the system may be 
expected to eventually suffer the phase transition into a 
Wigner crystal. 

To determine the relevant regime for cuprates one 
needs an estimate of the dimensionless coupling constant 
A. We find that 2ne 2 /ed w 2500K assuming e » 30 [Tt| 
and d = 12 A. in YBCO. The estimate of the temperature 




FIG. 3: The scaled B3| superfluid density K(T)/K(0) as 
a function of T/T c for various values of the dimensionless 
interaction A, and the data on underdoped YBCO [%^ . 



scale Td requires some assumptions on the relation be- 
tween the superfluid density of the bosons and the mea- 
sured superfluid density. For a finite interaction 



K(0) < — p , 

171 



(24) 



where the right-hand side is the undepleted superfluid 
density of the non-interacting system. In the effective 
gauge theories of the t- J model [6j , or of the fluctuating d- 
wave superconductor , density of bosons at low doping 
equals the density of holes, and thus po = x/a 2 , with a 
the lattice constant in the ab-plane. So 



T d > 



{a/df 



K(0). 



(25) 



Since (a/d) 2 as 0.1 in YBCO it is convenient to take 
the doping x — 0.1 as the reference point, so that > 
K(0) at this particular doping. Furthermore, the total 
superfluid density is related to the superfluid density of 
the bosons via the Ioffe-Larkin rule in Eq. 1, written 
more precisely as 



Kr o l(T)=K- 1 (T) + K- 1 (T) 



(26) 



where K qp {T) = K qp (0) - a 2 (2 ln(2)/ir)(v F /v A )T + 
0(T 2 ) is the contribution from nodal quasiparticles, with 
a as the corresponding Fermi liquid parameter \u\. Ex- 
panding to the first order in T one finds 



K tot (T) = K tot (0) 



irvA 



T, 



(27) 



where the Z = K(0)/K qp (0). Measurements of the su- 
perfluid density and thermal conductivity lead to an es- 
timate Za w 0.8 at x w 0.1 and therefore 



K(0) = ^# « 5^(0), 



(28) 
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assuming conservatively that a — 1. Estimating T c w 
65K at x = 0.1 in YBCO, the interpolation of known re- 
sults on the penetration depth yields 1/A 2 (0) w 50//xm 2 
|l9). which expressed in Kelvins [2(J leads to K to t{0) ~ 
400K. This leads to the estimate 



T d > 2000K. 



(29) 



A similar value of the lower bound is obtained by re- 
peating the exercise at optimal doping. The value of the 
dimcnsionless coupling A in YBCO is then 



A < 1.25. 



(30) 



There are at least two reasons to suspect that the value 
of the interaction parameter A may lie significantly below 
our estimated upper bound in the last equation. First, 
for A ~ 1 K(0) would be well below the non-interacting 
value of h 2 po/m, which would in turn yield a lower value 
of A. Second, our estimate is evidently very sensitive to 
the value of Z in Eq. 28; assuming a < 1, for example, 
would bring Z closer to unity and significantly increase 
the energy scale Td, and thus decrease the value of A. So 
the Eq. 30 should be understood as a very comfortable 
upper bound, with A most likely laying well below it. 

The Josephson coupling t may be related to the super- 
fluid density along the c-axis K c (0) as 



t 



K c {0) 
K(0)- 



(31) 



Assuming that the measured superfluid density in very 
underdoped cuprates is dominated by the bosonic com- 
ponent, the above ratio is 10 -4 , an d ap pears to become 
doping independent at low dopings | 21| . So, we estimate 



K(Q) 



t 



10 4 



(32) 



in the underdoped regime. 

As the critical temperature changes from T c = 92K at 
optimal doping to T c — 9K in the extremely underdoped 
regime in YBCO, the afr-plane superfluid density changes 
by roughly two orders of magnitude |22| . Assuming a 
constant Josephson coupling in this range leads to the 
left-hand inequality in Eq. 19, A ^> t/K(0), being com- 
fortably satisfied by our estimates, and the Coulomb in- 
teraction may be safely considered to be effectively short- 
ranged. Allowing the Josephson coupling to also decrease 
with undcrdoping, which may be closer to reality, only 
strengthens the above conclusion. The right-hand in- 
equality, however, would not be quite satisfied for Awl. 



As we argued, however, this is only an upper bound, and 
A may in fact be significantly smaller. To see the effect 
of the coupling strength on the form of the superfluid 
density, we plot in Fig. $K(T)/K(0) from Eq. 14 for 
various values of A |23|. together with the experimental 
points on YBCO with 9K < T c < 22K. The best 
fit is achieved for A = 10~ 3 , although it is clear that 
any value A < 0.1 would be almost equally good. This 
supports the recent proposal Q by one of us that the 
superfluid density in very underdoped cuprates is essen- 
tially the Bose condensate of the non-interacting layered 
bosonic system. The reader should also remember that 
whereas Eq. 14 yields the correct temperature depen- 
dence at low temperatures, it does not include the criti- 
cal fluctuations within the critical region of width ~ AT C 
near T c . These are known to modify the superfluid den- 
sity into K{T) ~ (T c - T)"-» , with v xy w 0.67, 0, and 
thus contribute to an additional rounding of the curve 
K(T) from its non-interacting form. The absence of any 
such discernible critical region in the data in Fig. 3 ad- 
ditionally supports our suggestion that A in YBCO may 
be rather small. 

If A ~ 1, the system crosses from a 2D, short-range 
regime, for K(0)/t <C 1, to a 3D long-range regime, for 
K(0)/t ~ 1, and the regimes II, III, and IV from Fig. 
1, well separated for weak coupling, now overlap signifi- 
cantly. The regime I, however, even in this case remains 
wide and distinct. In fact, it is well known that in 4 He, 
which is a strongly interacting Bose liquid, the variation 
of the superfluid density with temperature is well de- 
scribed by Landau's two-fluid model, except in the criti- 
cal region. It thus seems likely that Eq. 14 would remain 
qualitatively correct over most of the temperature region 
even if the bosonic system is not quite weakly interacting, 
as long as it is reasonably far from solidification. 

In conclusion, we discussed the nature of the collec- 
tive modes in a strongly anisotropic bosonic superfluid 
with Coulomb interaction between bosons. In particular, 
the influence of the anisotropic dispersion of these modes 
on the temperature dependence of the superfluid den- 
sity was analyzed. Depending on relative values of the 
anisotropy and the interaction four different regimes can 
be discerned for a weak interaction. A crude estimate of 
relevant parameters for cuprates shows that the Coulomb 
interaction in underdoped YBCO may be considered to 
effectively be short-ranged and at least marginally weak, 
and the system to be quasi-two-dimensional. 

We are grateful to S. Dodge, P. Lee, and S. Sondhi 
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NSERC of Canada. 
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